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Abstract 

The leading order evolution of parton orbital angular momenta is treated in the 
probabilistic manner. As a result, the splitting probability matrix, while coincid- 
ing with the recent explicit calculations, is expressed in the terms of the standard 
spin-dependent and spin-averaged kernels. The implied relations between spin and 
orbital angular momenta evolution kernels may be interpreted as resulting from the 
ambiguity in the separation of the angular momenta to the spin and orbital contri- 
butions and are likely to be valid for higher orders and twists. Consequently, the 
orbital and total angular momenta may be considered as decoupling from the spin 
ones, the latter being the only elements of nucleon structure accessible via the hard 
scattering of polarized particles. In particular, the partition of the total angular mo- 
menta of quarks and gluons should be analogous to that of their momentum at any 
Q 2 , generalizing the earlier finding for asymptotic limit by Ji, Tang and Hoodbhoy. 
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The orbital angular momenta of partons are known to be the necessary counterparts 
of their spins QCD evolution [|TJ. In particular, they are responsible for the cancellation of 
the large gluon helicity generated during the evolution, as was confirmed by the explicit 
calculation [Q, generalized recently to the higher moments of the orbital angular momenta 
distributions 0> @||- At the same time, it remains unclear, which general properties of the 
anomalous dimensions, if any, are relevant for such a cancellation. 

The more general question is the ambiguity of the splitting of the total angular momen- 
tum to the spin and orbital parts [§, ||, [7]]. In particular, the choice of energy momentum 
tensor suggested by Belinfante || allows to present the total angular momentum like the 
orbital one. One may ask, how such an ambiguity is compatible with the QCD evolution. 

In the present paper, the orbital angular momenta evolution is rederived using the 



kinetic interpretation [g] of the evolution equations || |K| [TT[ which appeared to be 
especially useful for the analysis of the spin dependence [0, [13| . As a result, the orbital 
angular momentum evolution is expressed entirely in the terms of the spin-dependent and 
spin-averaged kernels. While the conservation of total angular momentum is automatic, 
the extra relation between spin-dependent and spin-averaged kernels 

1 i rl 

dxxAP Gq (x) = - / dxxP Gq (x), (1) 
2, Jo 

guarantees, that the evolution of total angular momenta is equal to that of orbital an- 
gular momenta. This proves the compatibility of the QCD evolution with the Belinfante 
procedure. We discuss the related constraints for the NLO and power corrections to the 
angular momenta evolution. The emerging picture is a sort of decoupling of the spin 
and orbital/total angular momenta. While the first manifest themselves in the hard spin- 
dependent processes, the latter are the by-products of the momentum distribution and 
therefore related to unpolarized scattering. In particular, the similarity between evolution 
of parton momenta and angular momenta should result in the similarity of their partition 
at any Q 2 , 

The key element of the kinetic approach is the treatment of the virtual gluon contri- 
bution as a flow of the partons from the given point, while the terms associated to real 
gluons are interpreted as a flow to that point. As a result, the conservation of the op- 
erator like vector current or energy-momentum tensor acquire a form of particle number 
conservation. 

Say, the evolution equations for the singlet case, after writing || the virtual corrections 
to the splitting kernels like 



P(z)-6(l-z) f 1 P(t)dt, (2) 
Jo 



and performing the change of variable t = y/x in the last integral (which allows to interpret 
the partons decrease from virtual corrections as flow to the points y < x, similarly to that 
from real corrections), may be written as a Master equation of the Boltzmann type 

d ^ x ^ x ^ = ^J d y( w v(y ^ x )yfj(y) - w ii( x -* y) x fi( x ))- ( 3 ) 
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Here i,j are labels for q, G, and the transition rates per unit time t = InQ 2 



w l] {y^x) = ^P l A x6{y > X \ (4) 
Z7r y y 



The probabilistic interpretation of the orbital angular momentum evolution appeared 
already in the pioneering papers |], 0] . Say, the splitting of the quark to quark and gluon 
should be accompanied by the production of the orbital angular momentum, balancing 
the helicity of gluon (as the quark helicity is conserved). To calculate the net effect, one 
should subtract the probabilities of having the gluon with negative and positive helicities. 
Note that just the same combination (up to a sign) appears when one is calculating the 
gluon-quark spin- dependent splitting kernel, with the momentum fraction 1 — x: 



PjfW + P%{1 -x) = -AP G9 (1 - x), (5) 

where the upper inices are showing that the orbital angular momentum is produced by 
the spin one. The expression in the l.h.s appeared because the produced orbital angular 
momentum is carried by the quark with a momentum fraction x and gluon with a momen- 
tum fraction 1 — x. It is possible to find the ratio of the quark and gluon orbital angular 
momenta by means of the entirely probabilistic reasoning as well. Namely, suppose that 
the quark momentum before splitting has only z component, while the momenta of the 
final partons are in the x — z plane. The z components of the orbital angular momenta 
are 



L\ = P.rg; \F Z = -P x r° y , (6) 

where the x components of quarks and gluon momenta are equal, up to a sign, due to 
the momentum conservation, and the effective spatial non-locality of quark and gluon 
production r q ' G is introduced. The latter is leading also to the appearance of the x 
components of the orbital angular momenta: 

LI = ~Py y , l£ = -Pfrf. (7) 

Due to the conservation of the x component of the angular momentum L q x = —L^, 
and, consequently; 

r q pG 

v - * fQ\ 
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Note that the contribution of the x components of the parton helicities to the angular 
momentum conservation is resulting in the appearance of the additional terms of the order 
of Px/Pz m the r.h.s. of ® which should be neglected at leading order. Substituting fl|) 
to (|jD one get the partition rule 

L\ _ P? _ 1 -x 

±j z 1 z Jb 
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coinciding with the one found |§ by the explicit calculation. 

Although our derivation is a classical one and included the effective non-locality r, the 
latter did not enter the final result, and so due to the correspondence principle the above 
derivation, based on the conservation laws, is applicable in the actual quantum case as 
well. This would result in the substitution of the nonlocality r by the derivative operator 
in the momentum space. The simultaneous use of the analogs of (f| |7|) would mean the 
similarity in the action of such a derivative to the (non-forward) matrix elements of the 
T 0x and T 0z components of energy-momentum tensor, required by the Lorentz invariance. 

Applying the partition rule (|9|), one may deduce from (^) 

P q L q s (x) = (x-l)AP Gq (l-x), 

P£ q S (x) = (x-l)AP Gq (x). (10) 

Let us pass to the contribution to the orbital angular momentum due to the gluon helicity. 
The case of the quark-antiquark pair production is the most simple one, as their helicities 
are opposite in the chiral limit and all the gluon helicity is converted to the (anti) quark 
orbital momentum: 

P q L G s (x) = (l-x)P qG (x), (11) 

where the 1 — x factor is coming form the partition ratio(|9|), which does not depend on 
the nature of the splitted and produced partons. 

Consider now the case of gluon-gluon splitting, when the helicities of all the partici- 
pating gluons should be fixed: 

Pt£(x) = (1 - x)(P^ G + -(x) + P+ G -\x) - P+ G ++ (*)), (12) 

where the partition factor is appearing again, and the helicity of the gluon with the 
momentum fraction x is listed second, while the incoming gluon helicity is first, as usual. 
Actually, all the quantities in the r.h.s. were already calculated in the pioneering paper 
P~Q|] as an intermediate step of the computation of the spin-dependent kernel. However, 



it is instructive to use instead the general expressions: 

Pgg(x) = P^-{x) + P^ G - + (x) + P^{x), 

AP GG (x) = P^-{x) - P^ + (x) + P^ G ++ (x) (13) 
and the symmetry conditions 

P^-{x) = Pab + ^ ~ x), P^ + {x) = P+T (1 - x). (14) 
Combining these equations one easily get 

P^(x) = (x- l)(AP GG (x) + AP GG (1 -x)- P GG {x)). (15) 
Finally, the splitting matrix has the form: 

p LS (r) _ -f &P Gq {x) -2n f P qG (x) \ 

l j AP G9 (x) AP GG (x) + AP GG (x) - P GG (x) )> 
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where x = 1 — x and the obvious dependence on the flavour number is restored. 

Note that all the kernels in ( |16|) do not have a singular parts. It is a common feature 
in the kinetic approach, where the latter are treated separately via (0). It is natural 
therefore that (|i"6[ ) do not contain a singular terms at all. Recall, that such a terms 
are emerging due to the infrared divergence and are appearing in the diagonal kernels 
only. The resulting expression exactly coincides with the recent explicit calculation || £§] . 
However, the properties of the orbital angular momentum evolution can be seen now in a 
more general way. 

Namely, the evolution of the sum of the quark and gluon angular momenta, generated 
by the quark and gluon helicities, has a form: 

^-(L q + L G ) = f 1 dx(x-l)(AP Gq (l-x)+AP Gq (x))- [ 1 dxAJ:(x) + 
at Jo Jo 

[ l dx{x - l){AP GG {x) + AP GG (1 -x)- P GG {x) - 2n f P qG {x)) ■ [' dxAG(x). (17) 
JO Jo 

Changing x — > 1 — x in the second term of the quark coefficient and substituting 1 — x — > 
1/2 for the gluonic coefficient, which is possible because of the symmetry of the integrated 
function with respect to the interchange of x and 1 — x, one get: 

-{L q + L G ) = - f 1 dxAP Gq (x) ■ C dxAE(i) - 
at Jo Jo 

f 1 dx{AP GG {x) - xP GG (x) - 2n f xP qG (x)) ■ t dxAG(x). (18) 

JO JO 

The coefficients of AE and AG coincide, up to a sign, with those appearing in the evolution 
of gluon helicity (while quark helicity is conserved), guaranteeing the conservation of 
the total angular momentum. Note that the terms in the orbital angular momentum 
evolution, proportional to the unpolarized splitting kernels, are canceled by the terms in 
the gluon helicity evolution, coming from the virtual corrections, which are the same for 
the spin-dependent and spin-independent case ©. 

To proceed further, one should note that the diagonal evolution of the orbital angular 
momenta, generated by themselves, also does not include any new ingredients. As it was 
shown p| for the first moment and generalized recently by Hagler and Schafer || for the 
higher moments, this coincides with the evolution of unpolarized densities, with the extra 
power of x, P^ L (x) = xPij(x), so that the diagonal evolution of orbital angular momenta 

is precisely the same as that of the energy momentum tensors: PE L = P^ T ^ 2 \ where 
the moment is labeled in brackets. 

Combining all these facts, another interesting property of evolution may be found. 
Namely, the evolution of the total angular momenta J q = ^AS + L q , J G = AG + L G 
acquires the simple form: 




Jo 1 dx(x - l)P qq (x) 2n f Jq 1 dxxP qG (x) 
Jo dxxP Gq (x) -2n f $ dxxP qG (x 




and is governed by the same matrix P ii(1) = P TT ( 2 ). This result was first mentioned in 
T4, IT3 as a result of explicit calculation M and was adjusted to the proportionality of 



the relevant operators. 
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One should note that the equality 

p JJ(l) = pLL(l) (2Q) 

may be also considered as a quantitative expression of the Belinfante invariance, allowing 
for the representation of the total angular momentum in the orbital form. To understand 
this fact further, let us study, which properties of the kernels are implied by (|20|). It is 
easy to find the following relation between the matrices for the spin and orbital angular 
momenta evolution: 

pLS(l) , pSSQL) _ pLLQ), pJJ(l) _ p TT(2K 
ij ij ij v ij ij /' 

P? L (1) = 0. (21) 

Here P mean the evolution kernels, where distributions are weighted with their spin 
factors, which practically is resulting in the substitution P^q(x) = 1/2PM (x), pQg(x) = 
2P^(x), Pg g s = 2P q L q s , pQq(x) = 2P^(x), all other entries of P(x) and P(x) being equal. 
Note that these general relations imply the conservation of the total angular momentum 
as well, as it is now the consequence of the momentum conservation. 

The general reasons for such an equalities are the symmetry properties of the kernels 
with respect to the interchange of x and 1 — x, and the extra relation, required for the 
validity of the qq and Gq components of (f2ll). 

f 1 dxxAP Gq {x) = l - j l dxxP Gq (x), (22) 
Jo 2 Jo 

where 1/2 in the r.h.s. is just the quark spin appeared when transformation to P is 
performed. This equality may be easily checked explicitly. One may however ask, are 
there any general reasons for it. Such an explanation would be especially desirable, 
because of the relation to the Belinfante invariance, mentioned above. 

As soon as the origin of fl2~T|) is rather general, these relations should be satisfied for 
the higher orders of perturbation theory (where the whole probabilistic interpretation and 
the symmetry properties are, generally speaking, lost), as well as for the higher twistsQ 

Of course, this hypothesis should be checked by further explicit calculations. In par- 
ticular, higher loop calculations of Pss and Pls within dimensional regularization contain 
the extra ambiguity due to the choice of 75 matrix, and so it may either cancel identi- 
cally in their sum, or imply some particular factorization scheme, preserving Belinfante 
invariance, like it happen for Adler-Bardeen [ |TB| , |IT| relation or conformal invariance |2"H[ . 



In what follows, we shall suppose the validity of (pT|) and outline its possible physical 
consequences. 

The partition of the total angular momentum should always follow the partition of 
momentum: 

Jg(Q 2 ) = J2f Io l dxxq f {x,Q 2 ) 

J G (Q 2 ) Jo 1 dxxG( Xl Q 2 ) ' 1 ' 



1 In particular, the two-loop contribution of axial anomaly to is cancelled by the relevant term in 
Pq , which may be identified by calculation of the gluonic matrix elements of quark orbital momentum 
01- 
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As a result, at moderate Q 2 gluons should carry about half of the total angular momentum, 
as was anticipated by X. Ji |]l5f| , starting from the similar evolution and the results of sum 
rule calculation [103]. The Belinfante invariance is the new argument in favor of such a 
situation and one might look for the general non-perturbative proof of equality (23), based 
on the related arguments. 

Such a behavior would mean, that spin and orbital/total angular momenta are less 
dependent, when one might expect from the very beginning. The orbital and total an- 
gular momenta are completely defined by the momentum distribution, manifested in the 
scattering of the unpolarized particles and are in this sense decoupled from the spin ones. 
Recall, that total angular momentum conservation, being the main link between spin 
and orbital angular momenta, is actually guaranteed by momentum conservation and Be- 



linfante invariance ([H]). At the same time, the role of spin angular momenta in their 
canonical form is coming from their relation to OPE and manifestation in spin- dependent 
processes. 

In such a case, one should not expect any new information about nucleon spin structure 
from Deeply Virtual Compton Scattering |TJ] PL Instead, one may consider the relation of 
J and T as an extra constraint for the relevant non-forward parton distribution, namely: 



Y, I dxxE(x, A) = 0, (24) 
/ 



o 



which is just the term potentially making a difference between momentum and total 



angular momentum partition [14]. As the first moment of this quantity 

/ 

Jo 



dxE(x, A) = F 2 (0) (25) 



describes the anomalous magnetic moment of nucleon, the smallness of its isoscalar part 
is making such a behavior rather natural. Moreover, Belinfante invariance is providing 
the natural explanation of that fact (c.f. [LT7], [7]])! At the same time, model approaches 



22|j may provide the estimate of the effective quark orbital angular momentum as well as 
the higher moments of the total angular momenta, where ( PT| ) are not valid, and which 
are a natural objects for the models applications. Note also, that orbital momentum in 
field theory may essentially differ from that in quantum mechanics and in the potential 
models. 

Let us also briefly discuss the implications of another constraint in (|2~T]) , P$l = 0. This 
would mean that the orbital angular momentum is never produced by the spin one. At the 
same time, the possible gluon fusion at low x |^3| should give rise to such a production. 



One should conclude, that either this effect is oscillating, changing sign to give a zero in 
(pl|) after taking the first moment, or absent at all for some reasons. Note also, that this 
constraint is more general in a sense that it requires the similar evolution of L and J, 
while both may differ from that of T. 

In conclusion, we have found a representation for the leading order orbital angular 
momentum evolution, where its probabilistic nature is fully clear. The relevant splitting 



2 Although this process was originally suggested for this purpose, its subsequent studies pT|j , estab- 
lishing a new rapidly developing branch of PQCD, lost this original connection. 



6 



kernels are expressed via the usual spin-averaged and spin-dependent ones and are suitable 
for the numerical calculations. The implied constraints for the kernels are related to 
the ambiguity of separation of orbital and spin degrees of freedom and are expected 
to be valid for the higher order and power corrections. The resulting picture of the 
nucleon spin structure is rather simple. The orbital and (averaged over x) total momenta 
are decoupled and completely determined by the unpolarized scattering. The polarized 
scattering is entirely related to spin angular momenta, which should be studied in details 
in the forthcoming polarization experiments; the latter would be a final judgment for 
the proposed picture. The theoretical checks of the latter should include the study of 
orbital angular momentum within NLO approximations, higher twists corrections and 
non-perturbative approaches. 



I am indebted to A.V. Efremov, E. Leader, B. Pire, P.G. Ratcliffe, A. Schafer and J. 
Soffer for stimulating discussions. This research is partly performed in the framework of 
the Grants 96-02-17631 of Russian Foundation for Fundamental Research and Grant iV° 
93-1180 from INTAS. 
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